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Abstract 

To assess the durability of structures, heat and moisture transport need to be 
analyzed. To provide a reliable estimation of heat and moisture distribution 
in a certain structure, one needs to include all available information about 
the loading conditions and material parameters. Moreover, the informa- 
tion should be accompanied by a corresponding evaluation of its credibility. 
Here, the Bayesian inference is applied to combine different sources of infor- 
mation, so as to provide a more accurate estimation of heat and moisture 
fields The procedure is demonstrated on the probabilistic description of 
heterogeneous material where the uncertainties consist of a particular value 
of individual material characteristic and spatial fluctuations. As for the heat 
and moisture transfer, it is modelled in coupled setting j^. 

Keywords: uncertainty updating, Bayesian inference, heterogeneous 
materials, Karhunen-Loeve expansion, transport processes 



1. Introduction 

There are many important factors limiting the service life of buildings. An 
appropriate reliability analysis needs to take into account uncertainties in the 
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environmental conditions as well as in structural properties. Thanks to the 
growth of powerful computing resources and technology, recently developed 
procedures in the field of stochastic mechanics have become applicable to 
realistic engineering systems. 

The most common methods quantifying uncertainties are the first- and 
second-order reliability methods (FORM/SORM [ij) computing the proba- 
bility of failure related to limit states. Nevertheless, modern sophisticated 
and highly nonlinear models lead to non-Gaussian higher-order probability 
density functions of model parameters and its response. The higher moments 
become interesting quantities to be estimated by probabilistic analysis such 
as stochastic finite element methods (SFEM), see e.g. [H for a recent review. 
SFEM is a powerful tool in computational stochastic mechanics extending 
the classical deterministic finite element method (FEM) to the stochastic 
framework involving finite elements whose properties are random. 

This paper is focused on the modelling of uncertainties in material pro- 
perties and investigates the influence of such uncertainties on structural re- 
sponse. When dealing with homogeneous materials, one obtains a simplified 
scenario for SFEM, where uncertain material properties are described by ran- 
dom variables, which are assumed to be spatially constant. In the field of 
heterogeneous materials modelling, a widespread approach is multiscale mo- 
delling based on homogenization as presented e.g. in 0, 0, 0|. Nevertheless, 
homogenization theories have rigorous foundations for materials with well- 
defined geometry and components described by simple constitutive laws as in 
the case of regular masonry or composite materials with periodic microstruc- 
ture. However, these techniques cannot be efficiently applied to materials 
with random microstructure such as in case of quarry masonry, see Fig. [TJ 




Figure 1: Example of quarry masonry 

Another possibility is casting the description of heterogenous material 
within the probabilistic framework, where uncertain material properties in 
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time and/or space are represented by stochastic processes and fields. The 
resulting problem can be then solved by different strategies. The most fa- 
mous group involves spectral stochastic finite element methods (SSFEM) 
[gj]. These methods of uncertainty quantification are focused on propaga- 
tion of uncertainty in the system properties through the numerical model 
of the system in order to estimate probabilistic quantification of structural 
response. It assumes the knowledge of probabilistic formulation of uncertain 
system properties. Nevertheless, this uncertainty is usually very high before 
the structure is built, but after the construction, practical measurements 
can be performed. With these observations, the probabilistic models can be 
updated, to give a more accurate and reliable estimate. To this goal, the 
appropriate techniques from the field of inverse analysis should be employed. 

The following section presents a brief introduction into the inverse ana- 
lysis. Section [3] is devoted to Bayesian inference suitable for probabilistic 
estimation of model parameters from noisy and limited data. Section H] is 
focused on the probabilistic description of heterogeneous materials proper- 
ties where particular material parameters are not spatially constant. Section 
presents the application of Bayesian inference to model of coupled heat 
and moisture transfer in heterogeneous material and the obtained results are 
concluded in Section |6j 

2. Inverse problems 

In computational mechanics one tries to model a real system A, where 
system parameters q, a loading f and a system response u are related as 

A(u;q)=f. (1) 

The goal here is to obtain the response of the system for given parameters 
and loading conditions. In the field of inverse analysis, the goal is to find 
the values of system parameters q corresponding to given loading conditions 
f and experimental observations z. Therefore, one uses the numerical model 
of the system A and derives a so-called observation operator Y mapping the 
response u given parameters q and loading f to observed quantities z 

Y(q,f) = z. (2) 

The subject of this work is concerned with the description of heat and 
moisture conduction in structural materials and the system parameters q 
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are related here to the material behaviour. Material parameters are usually 
determined in the context of a chosen experimental setup, where the loading 
conditions are fixed, hence, the loading f is assumed to be constant in the 
following text. 

When simulating some real experiment, the model response is usually 
not equal to measured data because of experimental errors or imperfection 
of numerical model itself. It is often difficult to distinguish these sources of 
errors and they are described together in error vector e, so Eq. (j2J) becomes 

Y(q) + e = z. (3) 

The most common way of estimating material properties is based on 
fitting the response of numerical model to the results of real experiments, 
see e.g. (9[ for a recent review of parameters identification strategies. This 
usually leads to an optimization problem, where the difference between the 
model response and measured data is minimized by an appropriate optimiza- 
tion algorithm. Nevertheless, the formulation of the suitable error function 



is not always trivial [lOj, [ll|, [12|. The resulting function is often multimodal, 



non-smooth or non-differentiable and some robust optimization algorithm 



must be used (see e.g. [13| for an illustrative example of an application of 
evolutionary algorithm to solve such problem). 

There are only several works on including uncertainties into the curve 



fitting-based approach to parameters identification as e.g. in [14|. In general, 
these approaches omit related uncertainties in measurements as well as im- 
perfections of the numerical model and also the preliminary knowledge about 
the material parameters coming from their physical meaning (e.g. Young's 
modulus must be positive, Poisson's ratio must lay within the —1 and 0.5, 
etc). 

3. Bayesian updating of uncertainties 

Bayesian inference is the statistical inference in which the experimental 
observations are not used as the only source of information, but they are used 
to update the preliminary probabilistic description of system - the so-called 



prior information - to give the posterior distribution [15] . Recall that in re- 
alistic applications, observations are noisy, uncertain and limited in number 
relative to the dimension or complexity of the model space. Also, the model 
of a system may have limitations on its predictive value because of its im- 
precision, filtering or smoothing effects. Taking into account all pertinent 
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uncertainties, the process of material properties estimation cannot lead to a 
single 'optimal' parameter set, but one has to find a probability distribution 
of parameters that represents the knowledge about parameter values. The 
Bayesian setting for the inverse problems offers a rigorous foundation for in- 
ference from noisy data and uncertain forward models, a natural mechanism 
for incorporating prior information, and a quantitative assessment of uncer- 
tainty in the inferred results summarizing all available information about the 
unknown quantity In addition, unlike other techniques that aim to regu- 
larize the ill-posed inverse problem to achieve a point estimate, the Bayesian 
method treats the inverse problem as a well-posed problem in an expanded 
stochastic space. 

The Bayesian approaches to inverse problems have received much recent 
interest, since increasing performance of modern computers and clusters en- 
ables exhaustive Monte Carlo computations. Among recent applications one 



can cite applications in environmental modelling [16[, hydrology [17j or heat 



transfer [18[. We review this approach briefly below; for more extensive 
introductions, see [0]. 

The main principle of Bayesian inference is casting the inverse problem in 
the probabilistic setting, where material parameters q as well as observations 
z and also the response of forward operator Y(-) are considered as random 
variables or random fields. Therefore, we introduce the following notation. 
We consider a set J? of random elementary events u together with cr-algebra 
5? to which a real number in the interval [0, 1] may be assigned, the proba- 
bility of occurrence - mathematically a measure P. 

In the Bayesian setting, we assume three sources of information and un- 
certainties, which should be taken into account. The first one is our prior 
knowledge about the model/material parameters q(a>), which is represented 
by defining the prior density function p q (q). Prior models may embody sim- 
ple constraints on q, such as a range of feasible values, or may reflect more 
detailed knowledge about the parameters, such as correlations or smoothness. 

Other source of information comes from measurements, which are vio- 
lated by uncertain experimental errors e(u). Last uncertainty arises from 
imperfection of the numerical model included in the observation operator 
Y(-), when for example our description of the real system A does not include 
all important phenomena and therefore the forward operator response Y(-,u) 
can be assumed as uncertain. The probabilistic formulation of Eq. fl3]) now 
becomes 

z = Y(q(w),S) + e(0). (4) 
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If modelling uncertainties u cannot be neglected, they can be described 
by conditional probability density p(z|q) for predicted data z and given model 
parameters q. If these uncertainties can be neglected, only model parameters 
q(w) and observations z(u)) remain uncertain. In practise, it is sometimes 
difficult to distinguish the imperfection of the system description A from 
measurement error e. Hence modelling uncertainties u can be hidden in 
measuring error e(u). Finally, for noisy measurements we define the last 
probability density p z (z). 

To update our prior knowledge about model parameters we must include 
measurements with our theoretical knowledge. Bayesian update is based on 
the idea of Bayes' rule defined for probabilities. Definition of Bayes' rule for 
continuous distribution is, however, more problematic and hence [H, Chapter 
1.5] derived the posterior state of information 7r(q,z) as a conjunction of all 
information at hand 

7r(q,z) = KJ? q (q)ft(z)p(z|q) , (5) 

where k is a normalization constant. 

The posterior state of information defined in the space of model param- 
eters q is given by the marginal probability density 

7r q (q) = E Q [n(q, z)] = np q (q) / p(z|q)p z (z) ¥(du) = np q (q) L(q), (6) 

Jo 

where i? is a set of random elementary events u and measured data z enters 
through the likelihood function L(q), which gives a measure of how good a for- 
ward operator Y(q) is in explaining the data z. Here, is the expectation 
operator averaging over Q. 

To keep the presentation of different numerical aspects of particular meth- 
ods clear and transparent, we focus here on a quite common and simple case, 
where modelling-uncertainties are neglected and measurement errors are as- 
sumed to be Gaussian. Then the likelihood function takes the form 

m = *exp (-1 (Y(q) - z) T C7bs (Y(q) - z)) , (7) 

where C b s is a covariance among measurements z. 

The primary computational challenge is extracting information from the 
posterior density 7r q (q) jl|. Most estimates take the form of integrals over 
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the posterior, which may be computed with asymptotic methods, determin- 
istic methods, or sampling. The deterministic quadrature or cubature may 
be attractive alternatives to Monte Carlo simulation at low to moderate di- 
mensions, but Markov chain Monte Carlo (MCMC) [l9[ remains the most 
general and flexible method for complex and high-dimensional distributions. 



4. Uncertainty in properties of heterogeneous materials 

In modelling of heterogeneous material, some material parameters are 
not constants, but can be described as random fields. It means that the 
uncertainty in the particular material parameter q is modelled by defining 
q(x) for each x G Q as a random variable q{x) : Q — > 1R on a suitable 
probability space (i?, J5^,P) in some bounded admissible region Q C M. d . As 
a consequence, q : Q x Q — > IR is a random field and one may identify Q with 
the set of all possible values of q or with the space of all real- valued functions 
on Q. Alternatively, q(x, u) can be seen as a collection of real-valued random 
variables indexed by x e Q. 

Assuming the random field q(x,u) to be Gaussian, it is defined by its 
mean 

(i q (x)=E[q(x,u))]= [ q(x,u)¥(du) (8) 

J Q 

and its covariance 

C q (x,x') = E[(q(x,u) - fi q (x))(q(x',uj) - fi q (x'))} 

(q(x,co) - ^(x))(q(x\u) - N {x')) P(du;) . (9) 



Q 



Some non-Gaussian fields may be synthesized as - usually nonlinear - func- 



tions of Gaussian fields |20|, |21[. For instance, most of the material param- 
eters cannot be negative, hence, the lognormal random field can be more 
suitable for their description. Then the lognormal random field q(x,u) for 
each material parameter is obtained by nonlinear transformation of a stan- 



dard Gaussian random field q g (x,u) as, given in [21 



q(x,u) = exp(yu 9 + a g q g (x,u)) . (10) 

The statistical moments \i g and a g can be obtained from statistical moments 
fi q and a q given for lognormally distributed material property according to 
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following relations: 




<rj = ln( 1+ p) I , N = \n N -\a 2 g . (11) 



In a computational setting, the random field and the numerical model 
must be discretized. If the parameter field q(x) can be adequately represented 
on a finite collection of points {a?™ =1 } e M 2 , then we can write both the prior 
and posterior densities in terms of q = (q(xi), . . . , q(x n )), where qi = q(Xi) 
are random variables usually correlated among each other. The vector q, 
however, will probably be high- dimensional and it renders MCMC explo- 
ration of the posterior more challenging. Therefore, the Karhunen— Loeve 



expansion (KLE) can be applied to dimensionality reduction |22[ . 

KLE is an extremely useful tool for the concise representation of the 
stochastic processes. Based on the spectral decomposition of covariance func- 
tion C q (x,x') and the orthogonality of eigenfunctions ipi, the random field 
q(x, uj) can be written as 

oo 

?(^w) = /j 9 (x) + J^HW4 (12) 

8=0 

where £(o>) = (. . . , . . .) T is a set of uncorrelated random variables of 
zero mean and unit variance. The spatial KLE functions ipi(x) are the eigen- 
functions of the Fredholm integral equation with the covariance function as 
the integral kernel: 

/ C q (x,x')ipi(x)dx' = <^ijji{x) , (13) 
Jg 

where q are positive eigenvalues ordered in a descending order. 

Since the covariance is symmetric and positive definite, it can be expanded 
in the series 

oo 

C q {x 1 x') = Y,^A{*)A{x'). (14) 

i=l 

However, computing the eigenfunctions analytically is usually not feasi- 
ble. Therefore, one discretizes the covariance spatially according to chosen 
grid points (usually corresponding to a finite element mesh). The resulting 
covariance matrix C q is again symmetric and positive definite and Eq. (TT3l 
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becomes symmetric matrix eigenvalue problem, see [23J, where the eigen- 
functions ifti(x) are replaced by eigenvectors i/^- The eigenvalue problem is 
usually solved by a Krylov subspace method with a sparse matrix approxi- 



mation. For large eigenvalue problems, [24[ proposes the efficient low-rank 
and data sparse hierarchical matrix techniques. 

For practical implementation, the series (fT2j) and f lT^|) are truncated after 
M terms, yielding the approximations 



M 



q{u) « fi q + J2\^Uu)if>i, (15) 
i=i 

M 

c q « j^^-^- ( 16 ) 

Such spatial semi-discretization is optimal in the sense that the mean square 
error resulting from a truncation after the M-th term is minimized. 

5. Uncertainty updating in coupled heat and moisture tranfer 

This section is devoted to application of previously described techniques 
to uncertainty updating in coupled heat and moisture transport in heteroge- 
neous material with uncertain structure as quarry masonry. In particular, we 
employ the model proposed by Kiinzel described by the energy balance 
equation 

if ft = VT t AV ^] + h v V T [5 p VWp sat (6)}] (17) 
and the conservation of mass equation 

= VT [ D - V ^ + V T [<W mat (#)}] , (18) 

where 9 is the temperature, ip stands for the moisture and H, A, etc are de- 
scribed below. The transport coefficients defining the material behaviour are 
nonlinear functions of structural responses - the temperature and moisture 
fields - and material properties. We briefly recall their particular relations jij: 

• Thermal conductivity [Wm _1 K -1 ]: 

A^A fl + 6 "7< 6 -^ V (19) 
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Evaporation enthalpy of water [Jkg ]: 



273 15^ (°- 167 + 3 - 67 ' lcr4e ') 



h v = 2.5008 • 10 6 ( — — ) (20) 



Water vapour permeability [kgm x s x Pa 1 ]: 

h = 1.9446- io-' 2 , ( , + 273 . 15r . (21) 

Water vapour saturation pressure [Pa]: 

(17 08$ \ 
: « • (22) 
234.18 + 0/ V ; 



Liquid conduction coefficient [kgm s 



n = 3.8— • io ("-^/-d • £ ^ 

Wf (b — ip) 1 



(23) 



• Total enthalpy of building material [Jm 3 ]: 

# = p s c s ^ . (24) 

More detailed discussion about transport coefficients can be found in [2"!. lijj] . 
Some transport coefficients defined by Eqs. (Ti~9l - (|24|) depend on a subset 
of the material parameters listed in Tab. [TJ The approximation factor b 
appearing in Eqs. ([19]) and (|23|) can be determined from the relation: 

b = °- 8{W8 ° ~ ^ , (25) 

U>80 — 0.8l0/ 

where wso is the equilibrium water content at 0.8 [— ] relative humidity. There- 
fore, b is not considered as a material parameter, while uigo is another material 
property to be determined. Finally, there are 8 material parameters listed in 
Tab. [T]to be estimated by updating procedure. 

The presented table also contains the prior information about material 
parameters in terms of the mean values /i q and the standard deviations a q . 
Their particular values are chosen with regard to values corresponding to 
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Parameter 








o q 


Wf 


[kgm" 


i 


free water saturation 


200 


40 




[kgm" 


- 3 ] 


water content at 0.8 [— ] relative humidity 


100 


10 


Ao 


[Wm- 


-i K -i] 


thermal conductivity of dry material 


0.3 


0.1 


&tcs 


["] 




thermal conductivity supplement 


10 


2 


n 


["] 




water vapour diffusion resistance factor 


12 


5 


a 


[kgm" 


-2 s -0.5] 


water absorption coefficient 


0.6 


0.2 


c s 


[Jkg- 




specific heat capacity 


900 


100 


Ps 


[kgm" 


- 3 ] 


bulk density of building material 


1650 


50 



Table 1: Mean values and standard deviations of material parameters 



materials used in masonry 26[. Other prior information is that all these 
material parameters cannot be negative and hence, they could be considered 
as lognormally distributed. To describe the uncertainty about parameters of 
heterogeneous material, we choose the covariance kernel of a corresponding 
random field. Since the material properties change in the space because of 
changes in material components, we assume that the spatial fluctuations of 
all parameters are equal. It is probably not the best description of a real ma- 
terial. Nevertheless, we are convinced that the full spatial correlation among 
material properties is more realistic then full spatial independence. For the 
sake of simplicity, we do not study here the case of arbitrarily correlated pa- 
rameters. Therefore, we assume the same normalized exponential covariance 
kernel for all parameters, 



C(x,x') 



exp 





X ^ 




x 2 - x' 2 


(- 






1x2 



X 



X2), x — (x\, x' 2 ) 



(26) 

where 1 X1 and l X2 are covariance lengths. We assume also that the expert is 
certain about correlation lengths 1 X1 =0.1 [m] and l X2 = 0.04 [m], but he is 
not sure about a particular distribution of phases in material. In practise, 
the correlation lengths can be determined by the image analysis of a given 
material 27|, common size of bricks in masonry etc. 

Utilizing the covariance kernel fT26|) . we compute particular realizations 
of standard Gaussian random field based on Karhunen-Loeve expansion 



u 



(27) 
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Here, the eigenvectors describe the fluctuation of material property within 
the studied domain Q. Since the fluctuations reflect here the distribution of 
phases in material, the random variables & are also considered same for all 
material parameters. Given the prior mean fi q and standard deviation a q 
for each material parameter in Tab. (U corresponding statistical moments \i g 
and a g can be derived from Eq. (jlip . These can be then applied to Eq. ffTUl) 
in order to obtain a lognormal random field for each material parameter. 
Nevertheless, such procedure implies that the expert is also certain about 
mean value and relative amplitudes (given by the prior standard deviation) 
of respective random fields. That will unlikely happen in practice. To make 
the example more realistic, we include the uncertainty in the mean values of 
particular random fields by adding one random variable £^0 f° r each material 
property and extending the Eq. fflUj) into 

/ M \ 

q{u) = exp I fig + a g £ gfi + a g ^ >/s&( w )^i ) • (28) 

As a result, the random field corresponding to particular material property 
can be now shifted independently to each other. For a sake of simplicity 
and to keep the number of random variables in reasonable bounds, we keep 
the amplitudes of respective random fields still closely related. The total 
number of random variables L to be updated within Bayesian framework 
now becomes L = M + W, where M is the number of terms in the truncated 
KL expansion and W is the number of material properties. 

As an example, we consider two-dimensional rectangular domain dis- 
cretized by FE mesh into 80 nodes and 120 elements. Its geometry together 
with the specific loading conditions are shown in Fig. |2j The initial temper- 
ature is 9- m = 14 [°C] and moisture (fi n = 0.5 [-] in the whole domain. One 
side of the domain is submitted to exterior loading conditions 9 cxt = 5 [°C] 
and Lp cxt = 0.5 [-], while the opposite side is submitted to interior loading 
conditions 9- mt = 24 [°C] and (p int = 0.8 [-]. 

Exterior 
loading 
condition 

#ext = 5 [°C] 

<A=xt = 0-5 [- 

Figure 2: Experimental setup 
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In order to reduce the number of random variables to be updated within 
the Bayesian inference, one needs to choose the number of eigenmodes M as 
small as possible, but high enough for a satisfactory description of parameter 
fields. The attention should be paid to the error in description of parameter 
fields as well as to the related error in the model response. Fig. [3] presents 
a comparison of an arbitrary realization of thermal conductivity field Xo(x) 
computed using all 120 eigenmodes and its approximation Xq(x) computed 
using only first 7 eigenmodes. 



A [Wm-iR" 1 ] [Wm-'K" 1 ] 




(b) 

Figure 3: Thermal conductivity field computed using (a) all 120 eigenmodes and (b) only 
first 7 eigenmodes 

When comparing the corresponding heat and moisture fields, the negli- 
gence of higher eigenmodes in description of input parameter fields is reflected 
by fluctuations of response fields. These are, however, relatively small com- 
paring to absolute values of temperature and moisture (see Fig. [6]). For 
better understanding, Fig. H] shows only the fluctuations of response fields, 
since we have subtracted the response fields (9(x) and <p(x)) corresponding 
to homogeneous medium. (The particular response fields correspond to time 
* = 200 [h].) 

One can conclude that the employed Kiinzel's model has a smoothing 
effect, because negligence of higher eigenmodes induced relatively high er- 
ror in the approximation of input parameter fields, while its impact to the 
fluctuations of response fields is almost vanishing. 

In order to choose an appropriate number of eigenmodes, a relative point- 
wise error of input fields averaged over all 120 finite elements and over 100 
independent random realizations can be computed according to 

"loo ^120^ 55) • (29) 
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Figure 4: Fluctuations of response fields: (a) temperature field and (b) moisture field 
obtained for complete input parameter fields; (c) temperature field and (d) moisture field 
obtained for approximated input fields 



A similar error can be also computed in terms of response fields. These errors 
as a function of the number of eigenmodes M involved in the description 
of input fields computed for three different choices of correlation lengths 
are depicted in Fig. |5j It can be seen again that the error in description 
of input fields is decreasing slowly, while the error in the response fields 
descends much faster due to the smoothing effect of the numerical model. 
Owing to the results in Fig. [51 we decided to perform the Bayesian inference 
with an approximation of input fields using 7 KLE modes, which provides a 
satisfactory accuracy. 

According to the formulation of lognormal input fields given in Eq. (|2"H|) . 
we use one random variable for each eigenmode involved (it is M = 7) and 
one random variable for each material properties (it is W = 8) in order to 
enable their relative shift. Hence, we have L = 15 random variables to be 
updated within the Bayesian inference. 

Due to the lack of experimental data, we prepared a virtual experiment 
based on simulation including all eigenmodes. The values of temperature and 
moisture are measured in 14 points shown at Figs. |6] (a) and (c) and at three 
distinct times given in Figs. [6](b) and (d). Hence, the observations d consist 
of 84 values. They were then perturbed by Gaussian noise with standard 
deviation for temperature <jg = 0.2 [°C] and for moisture a v = 0.02 [— ] 



14 




Figure 5: Relative mean point- wise error (a) of the input thermal conductivity field and 
(b) of the overall responses induced by KLE approximation based on M eigenmodes 
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Figure 6: Virtual observations: (a) and (c) spatial arrangement of probes; (b) and (d) 
temporal organization of measurements 



and in that way we produced 100 virtual measurements. Based on them we 
calculated the observation covariance matrix C b s appearing in the likelihood 
function, which has the Gaussian form shown in the Eq. ((7j). 

The Bayesian update was performed using Metropolis-Hasting algorithm 
and 80, 000 samples were generated in order to sample the posterior density 
([5]) over variables £ = (£i . . -Cis)- Then, one realization of parameter fields 
was computed for each sample vector The mean computed over all pos- 
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terior fields of thermal conductivity E[A U)P osterior] depicted in Fig. [7J (a) can 
be compared with the reference field in Fig. [3] (a) utilized for preparation 
of virtual experiment. Since the mean is not a sufficient descriptor of the 



-A 




(a) (b) 

Figure 7: A posteriori fields of the thermal conductivity (a) mean over 80,000 samples 
and (b) cut along the axis x at y = 0.03 [m] of a subset of 200 posterior (light grey) and 
200 prior samples (dark grey) and the reference field (bold black) 

a posteriori distribution, Fig. [7J (b) shows a cut of a subset of 200 posterior 
samples (light grey lines). The cut is driven along the axis x at y = 0.03 [m]. 
These posterior samples can be compared with the same number of fields 
computed for samples drawn from the a priori distribution (dark grey lines) 
and the cut through the reference field of Ao (bold black line). One can see 
that the a posteriori samples much better encompass the reference field than 
the a priori ones. 

During the updating process, the model responses obtained for the a pos- 
teriori samples were stored. Fig. [8] shows the difference between the reference 
response fields and the mean computed over posterior response fields, both 
at time t = 200 [h]. Comparing the differences in Fig. [8] with size of fluctu- 
ations of corresponding fields in Figs. H] (a) and (b), one can conclude that 
the differences are relatively small. 

Finally, the a posteriori and the a priori distributions of the responses 
can be compared with the reference response at Fig. [9j The figure shows the 
evolution of the temperature (a) and the moisture (b) in time at the FE node 
no. 4. It is evident from the figures that the a posteriori samples (light grey) 
are distributed in the very vicinity of the reference response (bold black), 
while the dispersion of responses computed for 200 samples drawn from the 
a priori distribution is very large. 
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Figure 8: Difference between reference field and the mean over the a posteriori samples 
(a) for temperature field and (b) for moisture field 
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Figure 9: Temperature (a) and moisture (b) evolution in time at FE node 4: comparison 
of a priori samples (dark grey), a posteriori samples (light grey) and reference response 
(bold black) 



6. Conclusions 

The presented paper deals with the Bayesian updating of uncertainty in 
properties of heterogeneous materials. The process starts with the a priori 
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highly uncertain information given by an expert about material properties, 
which enters as input into the chosen material model. The model employed 
here is the Kunzel's model of coupled heat and moisture transfer, which 
assumes highly nonlinear relation between 8 material parameters and the 
structural response as given in Eqs. ffl~9]) - ff24"|) . Beside the prior uncertainty 
in the value of material characteristic itself, another uncertainty arises when 
describing heterogeneous material with spatial fluctuations of such property 
In order to describe the fluctuations, one random field is assigned to each 
material property. For a sake of simplicity, the fully correlated random fields 
are assumed, but one random variable is added into the description of each 
parameter so as to allow their mutual shift. To limit the number of ran- 
dom variables describing the material, the random fields are approximated 
by Karhunen-Loeve expansion with 7 eigenmodes. Finally, to grasp all men- 
tioned uncertainties, 15 random variables are considered. 

Beside the prior information, the virtual experiment is prepared to sub- 
stitute real experimental observations. The Metropolis-Hastings algorithm 
is then employed to sample the posterior distributions of random variables 
combining the prior knowledge and the information obtained from measure- 
ments. Figs. [7] and |9] present the results of the Bayesian inference verification. 
It is shown that even for highly nonlinear model, the updating process leads 
to much more precise prediction of the material properties as well as model 
responses. 

The drawback of the described procedure is the high computational cost. 
One simulation of the presented experiment takes 3.8 [s] at Intel Core Duo 
Processor T9600 with 4GB RAM. Hence, the whole sampling lasted almost 
90 [h]. Therefore, our future work will be focussed on the acceleration of sam- 



pling procedure via approximation of the structural response (see e.g. [28] ) or 



approximation of the posterior density 29] by polynomial chaos expansion. 



Acknowledgement 

This outcome has been achieved with the financial support of the Czech 
Science Foundation, project No. 103/08/1531 and No. 105/11/0411 and 
the Ministry of Education, Youth and Sports, project No. MSM6840770003. 
Many thanks belong also to Prof. H. G. Matthies, Ph.D. from TU Braun- 
schweig for a lot of fruitful discussions through the whole work on this paper. 



18 



References 

[1] A. Tarantola, Inverse Problem Theory and Methods for Model Param- 
eter Estimation, Society for Industrial and Applied Mathematics, 2005. 

[2] H. Kiinzel, K. Kiessl, Calculation of heat and moisture transfer in ex- 
posed building components, International Journal of Heat Mass Transfer 
40 (1997) 159-167. 

[3] O. Ditlevsen, H. O. Madsen, Structural Reliability Methods, John Wiley 
& Sons Ltd, Chichester, England, 1996. 

[4] G. Stefanou, The stochastic finite element method: Past, present and fu- 
ture, Computer Methods in Applied Mechanics and Engineering 198 (9- 
12) (2009) 1031-1051. 

[5] J. Sykora, J. Vorel, T. Krejci, M. Sejnoha, J. Sejnoha, Analysis of cou- 
pled heat and moisture transfer in masonry structures, Materials and 
Structures 42 (8) (2009) 1153-1167. 

[6] R. Valenta, M. Sejnoha, J. Zeman, Macroscopic constitutive law for 
mastic asphalt mixtures from multiscale modeling, International Journal 
for Multiscale Computational Engineering 8 (1) (2010) 131-149. 

[7] J. Vorel, M. Sejnoha, Evaluation of homogenized thermal conductivi- 
ties of imperfect carbon-carbon textile composites using the mori-tanaka 
method, Structural Engineering and Mechanics 33 (4) (2009) 429-446. 

[8] R. Ghanem, P. D. Spanos, Stochastic finite elements: A spectral ap- 
proach, second revised Edition, Dover Publications, Mineola, New York, 
2003. 

[9] A. Kucerova, Identification of nonlinear mechanical model parame- 
ters based on softcomputing methods, Ph.D. thesis, Ecole Normale 
Superieure de Cachan, Laboratoire de Mecanique et Technologie (2007). 

[10] G. Lubineau, A goal-oriented field measurement filtering technique for 
the identification of material model parameters, Computational Mechan- 
ics 44 (5) (2009) 591-603. 



19 



[11] A. Kucerova, Methodology of field measurements filtering to maximize 
the correlation with material parameters, in: WCSMO-8, Laboratorio 
Nacional de Engenharia Civil, Lisboa, Portugal, 2009, pp. on CD-ROM. 

[12] M. Kuraz, P. Mayer, M. Leps, D. Trpkosova, An adaptive time dis- 
cretization of the classical and the dual porosity model of Richards' 
equation, Journal of Computational and Applied Mathematics 233 (12) 
(2010) 3167-3177. 

[13] A. Kucerova, D. Brancherie, A. Ibrahimbegovic, J. Zeman, Z. Bittnar, 
Novel anisotropic continuum-discrete damage model capable of repre- 
senting localized failure of massive structures. Part II: identification 
from tests under heterogeneous stress field, Engineering Computations 
26 (1/2) (2009) 128-144. 

[14] D. Lehky, D. Novak, Probabilistic inverse analysis: Random material 
parameters of reinforced concrete frame, in: Ninth International Con- 
ference on Engineering Applications of Neural Networks, EAAN2005, 
Lille, France, 2005, pp. 147-154. 

[15] M. C. Kennedy, A. O'Hagan, Bayesian calibration of computer models, 
Journal of the Royal Statistical Society: Series B (Statistical Methodol- 
ogy) 63 (3) (2001) 425-464. 

[16] E. Yee, F.-S. Lien, A. Keats, R. D'Amours, Bayesian inversion of con- 
centration data: Source reconstruction in the adjoint representation of 
atmospheric diffusion, Journal of Wind Engineering and Industrial Aero- 
dynamics 96 (10-11) (2008) 1805-1816. 

[17] J. Fu, J. J. Gomez- Hernandez, Uncertainty assessment and data worth in 
groundwater flow and mass transport modeling using a blocking Markov 
chain Monte Carlo method, Journal of Hydrology 364 (3-4) (2009) 328- 
341. 

[18] S. Parthasarathy, C. Balaji, Estimation of parameters in multi-mode 
heat transfer problems using Bayesian inference effect of noise and a 
priori, International Journal of Heat and Mass Transfer 51 (9-10) (2008) 
2313-2334. 

[19] L. Tierney, Markov chains for exploring posterior distributions, Annals 
of Statistics 22 (4) (1994) 1701-1728. 



20 



[20] H. G. Matthies, Encyclopedia of Computational Mechanics, John Wi- 
ley &: Sons, Ltd., 2007, Ch. Uncertainty Quantification with Stochastic 
Finite Elements. 

[21] B. Rosic, H. G. Matthies, Computational approaches to inelastic media 
with uncertain parameters, Journal of the Serbian Society for Compu- 
tational Mechanics 2 (1) (2008) 28-43. 

[22] Y. Marzouk, H. Najm, Dimensionality reduction and polynomial chaos 
acceleration of Bayesian inference in inverse problems, Journal of Com- 
putational Physics 228 (6) (2009) 1862-1902. 

[23] H. Matthies, A. Keese, Galerkin methods for linear and nonlinear elliptic 
stochastic partial differential equations, Computer Methods in Applied 
Mechanics and Engineering 194 (12-16) (2005) 1295-1331. 

[24] B. N. Khoromskij, A. Litvinenko, Data sparse computation of the 
Karhunen-Loeve expansion, in: Proceedings of International Conference 
on Numerical Analysis and Applied Mathematics 2008, Vol. 1048, 2008, 
pp. 311-314. 

[25] R. Cerny, J. Madera, J. Kocf, E. Vejmelkova, Heat and moisture trans- 
port in porous materials involving cyclic wetting and drying, in: Com- 
putational Methods and Experimental Measurements XIV, Vol. 48 of 
WIT Transactions on Modelling and Simulation, 2009, pp. 3-12. 

[26] Z. Pavlfk, J. Mihulka, J. Zumar, R. Cerny, Experimental monitoring 
of moisture transfer across interfaces in brick masonry, in: Structural 
Faults and Repair, 2010. 

[27] M. Lombardo, J. Zeman, M. Sejnoha, G. Falsone, Stochastic modeling 
of chaotic masonry via mesostructural characterization, International 
Journal for Multiscale Computational Engineering 7 (2) (2009) 171-185. 

[28] A. Kucerova, H. G. Matthies, Uncertainty updating in the description of 
heterogeneous materials, Technische Mechanik 30 (1-3) (2010) 211-226. 

[29] B. Rosic, A. Litvinenko, O. Pajonk, H. G. Matthies, Direct Bayesian 
update of polynomial chaos representations, Journal of Computational 
Physics (0) (2011) Submitted for publication. 



21 



